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Abstract: A new method for analyzing arterial blood pressure is presented in this report. 
The technique is based on the scattering transform and consists in solving the spectral 
problem associated to a one-dimensional Schrodinger operator with a potential depending 
linearly upon the pressure. This potential is then expressed with the discrete spectrum 
which includes negative eigenvalues and corresponds to the interacting components of an N- 
soliton. The approach is similar to a nonlinear Fourier transform where the solitons play the 
role of sine and cosine components. The method provides new cardiovascular indices that 
seem to contain relevant physiological information. We first show how to use this approach 
to decompose the arterial blood pressure pulse into elementary waves and to reconstruct it 
or to separate its systolic and diastolic phases. Then we analyse the parameters computed 
from this technique in two physiological conditions, the head-up 60 degrees tilt test and 
the isometric handgrip test, widely used for studying short term cardiovascular control. 
Promising results are obtained. 
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Nouveaux indices cardio-vasculaires bases sur une 
analyse spectrale non lineaire de la forme des ondes de 

pression arterielle 

Resume : Une nouvelle methode d'analyse du signal de pression arterielle est proposee 
dans ce rapport. La technique est basee sur la transformation scattering et consiste a 
rcsoudrc un problcmc spectral associc a I'opcratcur dc Schrodinger unidimensionnel avec un 
potentiel qui depend lineairement de la pression. Le potentiel est reconstruit a I'aide du 
spectre discret uniquement qui comprend des valeurs propres negatives correspondant a N 
solitons en interaction. L'approche est analogue a une transformee de Fourier oil les solitons 
jouent le role des composantes sinus et cosinus. La methode introduit de nouveaux indices 
cardio-vasculaires qui semblent contenir des informations physiologiques interessantes. Dans 
un premier temps, nous appliquons la methode a la decomposition du signal de pression 
arterielle en ondes elementaires et a sa reconstruction ou a la separation de ses parties 
diastolique et systolique. Nous analysons par la suite les parametres calcules a partir de 
cette technique dans deux conditions physiologiques generalement utilisees pour I'etude du 
controle cardio-vasculaire a court-terme: test d'inclinaison a 60 degres et exercice isometrique 
(handgrip). Des resultats prometteurs ont ete obtenus. 

Mots-cles : Pression arterielle, transformation scattering, solitons, indices cardio-vasculaires 
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1 Introduction 

The analysis of mean values and beat-to-beat variability of cardiovascular (CV) time series 
has been widely used as a noninvasive approach to study the control of the autonomic ner- 
vous system (ANS) on the CV function [2], [23 • CV time series usually give information in 
the frequency and amplitude domains. The frequency (or period) is given by the RR interval 
(between two R peaks on the ECG ) or the pulse interval (PI) (between two systolic blood 
pressure peaks). The amplitude concerns systolic, diastolic and mean pressures (noted SBP, 
DBF and MBP respectively). Standard measures of these parameters are mean levels and 
global variability, spectral, temporal and time- frequency analysis [33] . 

Instead of the usual decomposition of the arterial blood pressure (ABP) waveform into 
a linear superposition of harmonic waves (sine and cosine functions) [ID] , [JB] , in this article 
we propose to use nonlinear superpositions of particular travelling waves, the N-solitons- 
solutions of the Korteweg-de Vries (KdV) equation where N travelling components are in- 
teracting. These N-solitons play the role here of the harmonic-waves solutions of the linear 
wave equation. The concept of soliton refers in fact to a solitary wave emerging unchanged 
in shape and speed from the collision with other solitary waves |38| . They fascinate scien- 
tists by their very interesting coherent-structure characteristics and are used in many fields 
to model natural phenomena. Solitons are solutions of nonlinear dispersive equations like 
the KdV equation arising in a variety of physical problems, for example to describe wave 
motion in shallow water canals [39] j [IH]- The use of solitons to describe the ABP was 
already introduced in [48 and '35] where a KdV equation and a Boussinesq equation were 
respectively proposed as a blood flow model. Recently, in |S], [21] a reduced model of the 
ABP cycle was introduced. The latter consists of a sum of a 2 or 3-soliton solution of a KdV 
equation, describing fast phenomena during the systolic phase and a 2-element windkessel 
model describing slow phenomena during the diastolic phase. We recall that the systolic 
phase corresponds to the contraction of the heart, driving blood out of the left ventricle 
while the diastolic phase corresponds to the period of relaxation of the heart. 

The decomposition of the ABP signal into a nonlinear superposition of solitons intro- 
duced in this article is based on an elegant mathematical transform: the scattering transform 
for a one-dimensional Schrodinger equation [5] , |12| , |14| . The main idea in our utilization of 
this transform consists in interpreting the pressure as a potential able to attract or repulse 
"fluid particles" or equivalently to transmit or reflect waves associated with them. This 
situation is modelled by a one-dimensional Schrodinger operator with a potential depending 
linearly upon the pressure wave [25] . The discrete levels of energy or speed of this system 
are given by the discrete spectrum of the Schrodinger operator. The associated eigenstates 
describe some coherent structures that are present in the pressure waveform and expressed 
as N solitons. 

The scattering-based signal analysis (SBSA) method introduces a new spectral descrip- 
tion of the pressure waveform leading to new cardiovascular indices. This study aims to 
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analyse these new parameters in two physiological conditions that are widely used for study- 
ing the short term control of the CV system: the head-up 60 degrees tilt test, in 15 healthy 
subjects, and the isometric handgrip exercise, in 13 healthy subjects. 

Among the SBSA parameters, some invariants of the scattering transform seem to be 
related to the contraction strength of the left ventricle and the stroke volume (SV). More- 
over, the beat-to-beat relation between the first eigenvalue and the heart period might be 
more reliable than the baroreflex slope for distinguishing two conditions. 

In the next section, we present the basis of the SBSA method. Section III compares real 
and reconstructed pressures using the SBSA and shows how we can separate the fast and 
slow pressure components related to the systolic and diastolic phases respectively. Section IV 
introduces the new cardiovascular indices computed using the SBSA technique and presents 
the results of the analysis in two physiological conditions: tilt and handgrip, including a 
discussion. Finally a conclusion summarizes the different results. 



2 A scattering-based signal analysis method 

In this section, we introduce a new signal analysis method based on the scattering theory. 
We start by briefly recalling the basis of the Direct and Inverse Scattering Transforms (DST 
& 1ST). Then, we present the main idea in the SBSA technique. For more details about 
DST and 1ST the reader can refer to the abundant literature and the references given. 

2.1 Scattering transform for a Schrodinger equation 

Let F be a given real function in the so-called Faddeev class L\{M.) [T]: 

/ + 00 
|l/(a;)|(H- |a;|)da; < oo}. (1) 
-OO 

We consider the one dimensional Schrodinger operator H{V) with a potential V\ 

H{V) : ^ ^ H{V)^ = -0 + V^. (2) 

The DST of V will be defined as a function of the solution of the spectral problem for H{V) 
where A and ip are respectively the eigenvalues and the associate eigenf unctions for some 
normalization: 

H{V)il: = (3) 

The spectrum of H{V) has two components: a continuous spectrum equal to (0, +oo) and 
a discrete spectrum with negative eigenvalues [TT], [T^], [2]. 
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For the positive eigenvalues denoted here A = fc^, there are eigenfunctions, called scatter- 
ing solutions of ^ consisting of linear combinations of exp (ikx) and exp {—ikx) as x —> ±00. 
Among them, we consider the Jost solutions from the left /; and from the right /r, normal- 
ized at ±00: 

(4) 
(5) 

(6) 
(7) 

(8) 

We recall that, for each fixed x e M the Jost solutions have analytic extensions in k to the 
upper- half complex plane [T] . 

The transmission coefficient T and the reflection coefficients Ri and Rr from the left and 
from the right respectively are defined through the relations: 

„ ,, , exp (ifca;) RAk) cxn i—ikx) ,„, 
fi{k, x) = ^1^^ ' + ' /^^^ '- + 0(1), X ^ (9) 

. , exp(-ifca:) Rr(k) exp (ikx) 

frik,x) ^ ^^^^^ + - + o{l), x^+00, (10) 

These coefficients satisfy: 

|T(fc)|2 + \Ri{k)\^ = |T(fc)|2 + |i?,(fc)|2 ^ 1. (11) 

On the other hand, for the negative eigenvalues of the discrete spectrum, equation Q 
admits solutions called bound states that belong to L^(R) in the x variable. When V belongs 
to the Faddeev class, the bound states solutions of ([3]) decay exponentially as x ±00 and 
their number N is finite [T]. Let us denote A„ = — with Ai < A2 < ... and ipn the N 
negative eigenvalues and L^-normalized bound states: 

/+00 
\Mx)\''dx = l, n = l,---,N. (12) 
-C30 

The eigenspaces being of dimension 1, the bound states and the Jost solutions are propor- 
tional: 

Ipnix) = CinfliiK„,x) = (-l)^""c™/r(»K„, x), (13) 

where c;„ and Cm are called the bound-state norming constants and are defined by: 

•-jn •- L / \fjiiK„,x)\'^dx]~^ , j = l,r. (14) 
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We now define the DST of V as tlie sets of scattering data from the left, Si{V) or from the 
right, SriV): 

«„, cj„, n = ,iV}, {jj}^{l,r}. (15) 

The potential can be uniquely reconstructed by using any one of these sets. The solution 
of this inverse problem, called 1ST, is the object of many studies concerned with specific 
classes of potentials [1], [7], [11], [15], [29]. Two transforms are then available, S and 
(in the sequel we choose j = r and drop the subscripts r and / for simplicity). 

In this study, we will use the special class of reflectionless potentials for which the left or 
right reflection coefficients are zero. Such potentials can be constructed as follows: let 
be the projector zeroing the i?-component of S{V), then o n^; o S{V) is reflectionless 
for any V in the Faddeev class. There are useful explicit representations of reflectionless 
potentials using only the discrete spectrum as in the following theorem [T3]: 

Theorem: If V is reflectionless for H{V)^ then: 

N 

V{x) ^ -A^Kni^l{x), XGM. (16) 

n=l 

V can be also written: 

^ a^(log(det(/ + A))) 
V(x) = -2 — , xeR, (17) 

where A is an x iV matrix: 

c c 

A{x) ^ [ — exp ((k„, + K„)x)], n,m=l,---,N. (18) 

Note that in (jl7p and (fTH|) , the potential is entirely defined with 27V parameters namely k„ 
and c„, 71 = 1, • • • , iV. 

A very close relation between a soliton solution of a KdV equation and a reflectionless 
potentiel of the Schrodinger operator was introduced in [Hj . In fact these potentials remain 
reflectionless when evolving in time and space according to a KdV equation. For t +oo, 
TV 1-solitons emerge, each one being characterized by a pair (k„, c„) such that gives the 
speed of the soliton and c„ its position. Therefore each component — 4k„7/;^ in the sum (|16p 
refers to a single soliton. 

2.2 A scattering based signal analysis method 

We now present how to use 1ST in a seemingly new method to analyse pulse-like signals of 
the Faddeev class. 
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The main idea in the SBSA approach is to interpret a positive signal y in the Faddeev 
class as a quantum well by changing the sign, and to tune the depth of this well with a 
positive parameter x i^i order to approximate y by a coherent state y^. For a deeper well 
the trapped energy will be higher and the approximation better, as we will prove. The 
estimate is then obtained by filtering out the nonlinear reflections: 

y^^ --S-^ oli^o S{-xy). (19) 
A convenient explicit formula is available, xVx being a reflect ionless potential: 

4 

^ n=l 

where —n\^n ^^-'^ V'x,"' — 1- • ■ • ■, are the negative eigenvalues and the associated L^- 
normalized eigenfunctions for H{—xy)- 

Then, we look for a value x for tJ^c parameter x such that the signal y is well approxi- 
mated by j/^. This is the decomposition of the signal y into the nonlinear superposition of 
solitons annomiced in j25j . 

It is well-known that the number of negative eigenvalues of H{~xy) is a nondecreasing 
function of x and there is an infinite unbounded sequence (xn) such that N^^ = N^^_^ + 1 
|25| , [31| . Determining the parameter x determines the number of negative eigenvalues and 
hence the number of solitons components required for a satisfying approximation of the 
signal y. Fig. [T] summarizes the SBSA technique. 



2.3 SBSA and the invariants of KdV 

The scattering transform has an infinite number of invariants which are related to the KdV 
conserved quantities [H], [52]. Let us denote these invariants Im{V), m = 0, 1, 2, • • • . They 
are of the form (we take — ^ as argument having (|19p in mind): 

i^i^v) = i~ir^'^-^£^ puv, )d^, (21) 

where Pm, m = 0, 1, 2, • • • are known polynomials in V and its successive derivatives with 
respect to a; G E [6J. 

A general formula relating /„(— V^), with — —xUi to the scattering data of H{V^) 
can be deduced; see for example [B], (TB], [52] : 

Im{xy)=T.4T'+^r^ {-k'r\n{\T^ik)\)dk, (22) 
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Figure 1: Signals analysis with the SBSA method 



m = 0, 1,2,---. 

We introduce the Riesz means of the negative eigenvalues A„ of H{V) such that A„ < 
A < 0: 

S'rA-y)= E ^^0- (23) 

A„<A 



Remark that So,\{V) is the number of eigenvalues of H{V) smaller than A. 
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For an A^Tj,-soliton, for instance —XVx ^^"^ previous subsection, the invariants only 
depend on the discrete spectrum and they are related to the Riesz means as follows: 

IniixVx) ^ S^fl{xVx)^ 7 = "^+^, m = 0, 1,2,--- (24) 
A "sum rule" is then verified by the invariants of xv a-nd xVx- 

Uxy) = Uxyx)+^^J J~k'r\n{\T^{k)\)dk, (25) 

m = 0, 1,2, • • •. 

In this article we are only interested in the two first invariants (m = and m = 1) 
corresponding to the conservation of mass and momentum for the KdV flows. Here it is 
sufficient to see them as invariants of the DST, in the same manner energy is invariant for 
the Fourier transform (Plancherel's theorem). We will show later in the application of the 
SBSA to the ABP that these two invariants are related to some important cardiovascular 
parameters. 

So, for m = 0, PoiVx, ■■■)^V^,we get with ^ and jSS]): 

ydx^ y^dx+ — / \n{\T^{k)\)dk. (26) 

) J —OO J~oo 

For m = 1, Pi(Fx, ■■■) = V^, we have with jUl) and jH]): 

/ + C>C o /' + OC 

vldx k^\n{\T^{k)\)dk. (27) 

-oo J — oo 

Equation (|27p is known as the Buslaev-Faddeed-Zakharov trace formula. 

Proposition: Let y : R ^ R be a continuous non-negative function with a compact 
support, then we have the convergence of the estimates of the first two invariants: 

lim Iraiyx) = I-m{y), "1 = 0,1. (28) 

X— »+oo 

Proof: We can apply the results on the Lieb-Thirring semiclassical limit of the Riesz means 

i, m- 

S^M ^ ^ci I ^ > 1 (29) 

X-++00 ' 2 

where L'^ is the so-called Lieb-Thirring constant given by: 



i^;'-(4-)-^^^. (30) 

r(7+2) 
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We notice that for 7 = -: 

lim ^' = I v{x)dx, Lf = -. (31) 

So, we deduce the convergence of the first invariant estimate. 
3 

For 7 = — we have an analog of the Plancherel identity for the Fourier transform: 

lim ^i4^ = Lf I y{xfdx^ Lf = ^ (32) 
Therefore, we get the convergence of the second invariant. 

3 Application of the SBSA to the ABP waves 
3.1 ABP reconstruction 

In the previous subsection, we presented a new signal analysis method based on the scatter- 
ing transform. Now, we propose to use this method for ABP analysis. For convenience we 
replace the space variable x by the time variable t. 

We note the ABP signal P{t) and the estimated pressure with the SBSA technique P{t) 
such that: 

4 

^W = -E«^x,«<„W, (33) 

n=l 

where — k^,„ and V'x."' " — I'''' '^x ^"^^ negative eigenvalues and associated 

L^— normalized eigenf unctions of H{—xP)- We recall that each component 4k^^„?/'^ „ in 
refers to a single soliton [H], 

In Fig. [2] and Fig. [3l measured and reconstructed pressures at the aorta and at the finger 
levels are presented respectively. The aortic pressure was measured using a catheter while a 
Finapres was used to measure the pressure at finger. Only 5 to 10 components are sufficient 
for a good reconstruction of the ABP waveform. 



3.2 Separation of the systolic and diastolic phases 

Following the work done in [9l , [24] , we propose here using the SBSA technique to separate 
the pressure into fast and slow parts corresponding respectively to the systolic and diastolic 
phases. Indeed a reduced model of ABP has been proposed in [S], [21]. The latter consists of 
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Figure 2: Measured and reconstructed pressure at the aorta 
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Figure 3: Measured and reconstructed pressure at the finger 



a sum of two terms: a 2 or 3-sohton solution of a KdV equation describing fast phenomena 
which predominate during the systohc phase and a 2-element windkessel model describing 
slow phenomena during the diastolic phase. As noticed in the previous section, the SBSA 
technique decomposes the ABP signal into a sum of solitons, each one characterized by its 
velocity given by the discrete eigenvalues —n\^n- So the largest k^,„, n = 1, • • • ,Ns describe 
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fast phenomena while the smallest ones describe slow phenomena. Referring to [S], [Mj, 
we take A^s = 2 or 3. We note Ps and Pd the estimated systolic and diastolic pressures 
respectively such that: 

P^^^) = ^Jl^X,n^ln. = ^ E ^X,n^ln- (34) 

n=l ^ n=Ns + l 

We compute the first two invariants of these partial pressures with the Riesz means for the 
chosen cut-off speed A: 

4 16 

INV,{X)^-S.^^{xP). INV2{\)^-^Ss^^{xP), (35) 

We can now define the proposed invariants for the whole beat and for the systolic and 
diastolic phases (INVj, INVSj, INVDj, j = 1,2 respectively): 

INVj = INVj{0), 

INVS, = INVjiK), j = 1, 2, (36) 

INVDj = INVj{0) - INVSj{\s), 

where — A^^2 or A^^a. 

In Fig. [4] and Fig. El we represent the measured pressure and the estimated systolic and 
diastolic parts respectively. We remark that Ps and Pd are respectively localized during the 
systole and the diastole, as expected. 



4 New spectral cardiovascular indices 
4.1 The SBSA parameters and the ABP 

The SBSA technique provides a new description of the ABP signal using the DST. As seen 
in the previous section, the reconstruction of the signal by 1ST from its spectral data gives 
good results. Instead of reconstructing the original signal, it is possible to modify the spec- 
trum leading to some kind of filtering. This is illustrated by the separation of the systolic 
and diastolic phases. Moreover, the SBSA method introduces new parameters. The first 
two global invariants INVi and INV2 are respectively, by definition, the usual mean blood 
pressure (MBP) and the less usual, but easy to compute directly, integral of the square of 
the pressure. The first systolic invariant INVSi is a new index: we think that it can be 
correlated to SV. Remark that if INV2 is easy to compute directly, the "fast part" of this 
integral, the second systolic invariant INVS2 is a new less obvious index and might con- 
tain information on ventricular contractility. SV and contractility are in fact parameters of 
great interest that are difficult to measure routinely, as they require invasive or sophisticated 
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Figure 4: Pg and fast systolic phenomena 




Figure 5: Pa and slow diastolic phenomena 



techniques. For instance SV can be estimated by invasive nuclear ventriculography (47j . 2D 
echocardiography [23] , radionuclide monitoring j43] , impedance cardiography [8] , [10] . Only 
one evaluation of SV from ABP has been proposed [TH] • Ventricular contractility is assessed 
by the mean of the tissue doppler echocardiography [T7] , [H] . 
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On the other hand, the eigenvalues computed with the SBSA technique are strongly de- 
pendent upon the ABP waveform. Indeed, when the pressure waveform changes, the optimal 
value of the parameter x and the eigenvalues also change. This fact is illustrated in Fig. 
[51 Therefore, the eigenvalues could be used to assess the baroreflex sensitivity (BRS) in a 
certain way. In fact, the BRS expresses the variation of the heart beat interval in response to 
each arterial pressure variation. The BRS concept was first based on drug-induced responses 
of ABP and heart period [27]. Then various time [27], [Si!, [SS] and spectral [S], [21], [5i] . 
domain methods [27j were compared. 

In this section we analyse the SBSA parameters in two situations, devoted to the assess- 
ment of the ANS control of the CV system. We restrict the study to the modulus of the 
first two negative eigenvalues lA^^i] = |Ai|, |A;^,2| = 1^21 and the first two invariants (global, 
systolic and diastolic). We include in the analysis some classical parameters which are PI, 
SBP and DBP. 



4.2 Head-up 60 degrees tilt-test 

The head-up tilt test is mainly used for vasovagal syncopes diagnosis, characterized by an 
autonomic dysfunction. It consists in the orthostatic transition from the supine to the 
standing positions. This leads to a redistribution of the venous blood volume, from the 
intrathoracic region towards the venous volume in the leg and lower abdominal veins. This 
leads to a decrease in SV and PI and an increase in SBP [TD] , [H] , [H] , [33] . 

A group of 15 healthy subjects under 0.25Hz paced breathing, already studied 3], was 
considered. The table was rotated to an upright position at 60 degrees. The continuous 
ABP was measured at the finger using a Finapres device [45j . The two positions, supine 
and standing, were compared using the Wilcoxon non parametric paired test. 

Fig. [7| shows the time series of PI, SBP, INVSi and |Ai| in the supine and standing 
positions. Mean levels of the ABP parameters are presented in Table \T\ We notice that 
significant differences between the supine and standing positions appear for PI, DBP, |Ai|, 
IA2I and INVSi while for the other parameters the differences are not significant. So after 
the tilt test, SBP and MBP recover their prior values (preceding the tilt phase) whereas 
INVSi is decreased. If INVSi is related to SV then this means that SV remains decreased 
after the tilt test and this can be explained by the decrease in the venous return. 

The ABP analysis during the transition from the supine to the standing positions was 
possible for only nine subjects. For these, we assessed the linear relation between PI of the 
beat n + 1 and |Ai| computed for the beat n. In Fig. [S] we note that as PI decreases, 
I All increases and that the linear correlation between PI and |Ai| is stronger around the 
transition than in supine position. Table [2] shows that the correlation coefficient (i?^) and 
the slope, compared by a one-way repeated measures analysis of variance, were significantly 
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Figure 6: (a) First eigenvalue time series of a healthy subjects under 0.25iJz paced-breathing 
during 29 beats. The signal exhibits a beat-to-beat variability that is twice as fast as the 
breathing rate, (b) A zoom on a few beats shows morphological changes in the blood pressure 
signal. 



stronger around the transition than in supine or standing positions. Moreover, a comparison 
with usual BRS indices, SBP and pulse pressure (PP), shows that |Ai| has the strongest 
correlation with PI (Table [31). It is not surprising that |Ai| is more informative about the 
relation between the heart period and the ABP because it reflects the arterial waveform and 
not only one (SBP) or two (PP) samples of this waveform. 

4.3 Isometric handgrip exercise 

The isometric handgrip exercise is mainly used for the evaluation of a non appropriate ANS 
behavior, that mimics actual conditions of professional or domestic exposure, with arterial 
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Figure 7: Time series of ABP parameters in a healthy subject under 0.25Hz paced breathing, 
in supine (a) and standing (b) positions after a 60 degrees tilt test. PI and INVSi are 
reduced whereas |Ai| is increased in the standing position, more than two minutes after the 
tilt test. SBP has the same level in both positions. 



hypertension. Indeed, the ANS acts in the same way as in a dynamic exercise usually char- 
acterized by an increase in muscle oxygen needs. The isometric exercise is a form of exercise 
involving the static contraction of a muscle without any visible movement in the angle of 
the joints [13j . 

A group of 13 healthy subjects was considered. The continuous ABP was measured at 
the finger with a Finapres device [45j . The two conditions: at rest and during the handgrip 
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Table 1: ABP parameters during tilt protocol in 15 healthy subjects 



ABP parameters 


supine 


standing 


probability 


Direct 








PI 


918 ± 33 


778 ± 21 




SBP 


121 ±4 


127 ±3 


NS 


DBP 


68 ±4 


80 ±3 




Eigenvalues 








First 


1515 ± 131 


1965 ± 137 




Second 


1205 ± 86 


1612 ±92 




First invariants 








Global 


77 ±6 


71 ±4 


NS 


Systolic 


21 ± 1 


19 ± 1 


** 


Diastolic 


56 ±4 


51 ±3 


NS 


Second invariants 








Global 


7041 ± 975 


6833 ±711 


NS 


Systolic 


2705 ± 328 


2566 ± 247 


NS 


Diastolic 


4336 ± 648 


4276 ± 465 


NS 



Data are expressed as means and SEM; **: P < .01; ***: P < .001 



Table 2: Beat-to-beat BRS during tilt protocol 





supine 


tilt 


standing 


probability 




.335 ± .102 


.611 ± .074 


.421 ± .079 




slope 


-.105 ± .03 


-.177 ± .03 


-.106 ± .015 





and slope of the linear regression between |Ai(n)| and PI(n + 1), over about 60 beats, in 9 healthy 
subjects. Data are expressed as means and SEM; ** p < .01; *** p < .001, at repeated measures analysis of 
variance. Slope and are the strongest during tilt. 



Table 3: Three indices of beat-to-beat BRS during tilt 





tilt 


|Ai| 


.611 ± .074 


SBP 


.351 ± .091 


PP 


.415 ± .069 



Mean of R^ of the linear regression between ABP parameters {|Ai|, SBP, PP) and PI, over about 60 beats, 
in 9 healthy subjects of the tilt test. Data are expressed as means and SEM. |Ai| is the most strongly 
correlated with PI. 

test, were compared using the Wilcoxon non parametric paired test. 

In Fig. [HI the time series of PI, SBP, INVS2 and |Ai| are presented at rest and during 
the handgrip test. Table|3]illustrates the mean levels of the ABP parameters. We notice that 
while all the first invariants do not change significantly, the second invariants are sensitively 
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Figure 8: Beat-to-beat BRS, represented as the relation between |Ai(n)l and PI{n + 1) 
Slope and correlation are stronger during tilt. 



increased during handgrip. 

The voluntary central command, involved in the handgrip, synchronously activates the 
motor and CV systems, leading first to an increase in the heart rate, followed by an increase 
in the ABP [31], [SU]. So, as expected, PI decreases and SBP increases (Table H}. Moreover, 
INVSi might inform us that SV does not change while INVS2 might inform us about the 
increasing contractility, as if the heart tries to eject the same quantity of blood during a 
smaller period. Such a result evokes the treppe effect (or frequency- force relation) , where an 
increase in heart rate indirectly induces an increase in contractility. 

As in the case of the tilt test, we study the relationship between PI and |Ai|. Fig. [TU] 
illustrates the relation between PI of the beat n + 1 and | Ai | computed for the beat n at rest 
and during the handgrip test. The strong linear correlation between PI and |Ai| is the same 
at rest and during the handgrip, but the slope is significantly lower during the handgrip 
(Table [K]). Moreover, a comparison with usual BRS indices, SBP and PP, shows that |Ai| 
has the strongest correlation with PI (Table EJ. This result, as obvious as in the case of the 
tilt test leads us to consider |Ai| as a promising index that can be used to study the relation 
between the heart period and ABP. 
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Figure 9: Time series of ABP parameters in a healthy subject under spontaneous breathing, 
at rest (a) and during the handgrip test (b). PI is reduced whereas SBP, INVS2 and |Ai| 
are strongly increased during handgrip. 



Table 4: ABP parameters during handgrip protocol in 13 healthy subjects 



ABP parameters 


rest 


handgrip 


probability 


Direct 








PI 


953 ± 59 


748 ± 42 


*** 


SBP 


148 ±5 


183 ±5 


*** 


DBF 


83 ±3 


104 ±4 


*** 


Eigenvalues 








First 


1439 ± 160 


2325 ± 249 


** 


Second 


1186 ± 125 


1872 ± 180 


*** 


First invariants 








Global 


96 ±6 


95 ±6 


NS 


Systolic 


26 ± 1 


26 ± 1 


NS 


Diastolic 


70 ±4 


69 ±4 


NS 


Second invariants 








Global 


10239 ± 866 


12886 ± 1192 


** 


Systolic 


3956 ± 344 


5013 ± 470 


** 


Diastolic 


6283 ± 526 


7872 ± 727 


** 



Data are expressed as means and SEM; ** p < .01; *** p < .001. 
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Table 5: Bcat-to-boat BRS durinp; handgrip protocol 





rest 


handgrip 


probability 




.647 ± .054 


.609 ± .050 


NS 


slope 


-.302 ±.071 


-.132 ±.024 


** 



and slope of the linear regression between |Ai(n)| and PI(n + 1), over about 40 beats, in 13 healthy 
subjects. Data are expressed as means and SEM; ** p < .01; *** p < .001, at paired test. Slope is reduced 
during handgrip whereas R? is not significantly diflferent. 



Table 6: Three indices of beat-to-beat BRS during handgrip protocol 





rest 


handgrip 


|Ail 


.647 ± .054 


.609 ± .050 


SBP 


.166 ±.055 


.261 ± .060 


PP 


.375 ± .078 


.284 ± .060 



Mean i?^ of the linear regression between ABP parameters (|Ai|, SBP, PP) and PI, over about 40 beats, 
in the 13 healthy subjects of the handgrip test. Data are expressed as means and SEM. is the most 
strongly correlated with PI. 
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Figure 10: Beat-to-beat BRS, represented as the relation between |Ai(n)| and PI{n + 1). 
Slope is lower during handgrip. 
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5 Conclusion 

This article deals with a new ABP analysis method based on the scattering theory. This 
SBSA method can be thought of as a nonlinear Fourier analysis for pulse-like signals. It 
allows analysis and precise reconstruction as is shown by the very good agreement between 
real and estimated pressures. We have also presented an application to a filtering problem 
consisting in separating the systolic and diastolic phases. Then, we have introduced new 
cardiovascular indices computed with the SBSA method. These parameters include the first 
two systolic invariants and we think that they might give information on the variation of 
the stroke volume and the ventricular contractility, that are difficult to measure routinely. 
Another interesting parameter is the first eigenvalue which seems to reflect the BRS in a 
certain way. The results obtained from the analysis of two widely used physiological condi- 
tions are promising and we are now working on the validation of the advanced hypotheses. 
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